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摘要 : 【目的 ] 利 用 核糖 体 DNA 联合 序列 探讨 天 牛 总 科 高 阶 元 分 子 系统 发 育 。【 方 法 】 本 研究 采用 分 子 标记 技术 ， 
分 析 测 定 了 63 种 天 牛 核糖 体 28S rDNA D2 和 D3 区 以 及 18S rDNA VA 和 V7 区 的 DNA 序列 ,并 采用 邻接 法 、 最 大 
似 然 法 和 贝 叶 斯 推论 法 分 别 构建 了 天 牛 总 科 2 科 6 亚 科 63 种 的 分 子 进 化 系统 。[【 结果 】 序 列 联 合 比 对 分 析 , 最 终 
得 到 1 404 bp 的 联合 数据 组 ,其 中 可 变 位 点 446 个 (32.0% ) ,保守 位 点 958(68.0% ) ,转换 / 丰 换 的 平均 值 (R 值 ) 为 
1.73, 288 rDNA 和 18S rDNA 以 及 联合 序列 的 饱和 度 分 析 显 示 碱 基 突 变 未 达到 饱和 ,说 明 这 些 序列 适合 于 分 子 进 
化 树 的 构建 。 利 用 不 同系 统 发 育 重 建 方法 得 到 进化 树 具 有 相似 拓扑 结构 ,结果 支持 沟 胜 天 牛 亚 科 、 花 天 牛 亚 科 和 
天 牛 亚 科 为 单 系 群 ,这 与 形态 学 分 类 结果 相似 ; 狭 胸 天 牛 独立 成 为 亚 科 得 到 了 支持 。【 结论 ] 利 用 28S rDNA D2 和 
D3 区 以 及 18S rDNA V4 和 V7 区 联合 序列 成 功 构 建 出 了 天 牛 总 科 高 阶 元 的 系统 发 育 树 。 人 研究 表明 联合 序列 分 析 
是 探讨 天 牛 高 阶 元 分 类 的 有 效 的 方法 。 
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Abstract: [ Aim] To explore the molecular phylogeny of higher taxa of the superfamily Cerambycoidea by 
using the combined analysis of ribosomal DNA sequences. [Methods] Ribosomal 28S rDNA ( D2 and D3 
region) and 18S rDNA ( V4 and V7 region) sequences from total 63 species representing six subfamilies 
and two families of Cerambycoidea were chosen as the molecular markers and further sequenced. 
Phylogenetic trees for the two families and six subfamilies of Cerambycoidea were constructed using 
Neighbor-Jointing ( NJ) , Maximum-Likelihood ( ML) and Bayesian-Inference ( BI) methods. [Results] 
Sequence analysis results indicated that 1 404 sites after alignment consisted of 446 variable sites 
(32.096 ) and 958 conserved sites (68. 0% ). The transition versus transversion ratio ( R-value) was 
1. 73. Saturation plot indicated that combined analysis had better linear relation than the single gene 
analysis, suggesting that the base mutations have not reached saturation. Phylogenetic trees based on NJ, 
ML and BI methods showed similar topology, suggesting that Lamiinae, Cerambycinae and Lepturinae 
were all monophyletic group, which was consistent with the results of the traditional morphology. Our 
results supported the subfamilial placement of Philinae. [Conclusion] We successfully reconstructed the 
phylogenetic trees of higher taxa of Cerambycoidea based on 288 rDNA ( D2 and D3 region) and 18S 
rDNA ( V4 and V7 region) sequences. This study confirms that the combined analysis of sequences is an 


effective method for studying the higher taxonomy of longicorn beetles. 
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K 4# P} ( Cerambycidae ) 3E Jj F 88 38 H 
( Coleoptera ) 天 牛 总 科 ( Cerambycoidea ) ,全 世界 已 记 
录 45 000 多 种 ,是 鞘翅 目 最 大 的 科 之 一 ( 华 立 中 等 ， 
2009) 。 天 和牛 总 科 的 分 类 经 历 了 广泛 的 修订 与 完 
善 ,已 形成 较为 完善 的 分 类 系统 (将 书 杭 等 ,1985 ; 
林 美 英 ,2011 ) Švácha 等 (1986 ,1997 ) 将 天 和 牛 科 分 
类 地 位 上 升 为 天 牛 总 科 , 并 提出 了 包括 以 下 各 科 的 
分 类 系统 :着 天 牛 科 、 盾 天 牛 科 、 上 暗 天 牛 科 和 天 牛 科 
(包括 异 天 牛 亚 科 、 饮 天 牛 亚 科 、 幽 天 牛 亚 科 、 椎 天 
牛 亚 科 、 花 天 牛 亚 科 、 膜 花 天 牛 亚 科 、 锯 花 天 牛 亚 科 、 
AIF MVRHRITI EE ACIE ERE) 。 但 是 天 牛 总 科 各 亚 科 
之 间 的 演化 和 归属 一 直 存 在 争议 ( 董 抗 震 和 杨 星 
FF, 2003) :Nakamura(1981 ) 根 据 天 和牛 幼虫 的 形态 学 
和 生物 学 人 研究 均 表 明 最 原始 的 为 寞 天 牛 亚 科 , 沟 用 
天 牛 亚 科 为 最 进化 的 ,其 分 化 顺序 为 异 天 牛 亚 科 、 饥 
天 牛 亚 科 一 瘦 天 牛 亚 科 一 花 天 牛 亚 科 一 椎 天 牛 亚 科 
一 幽 天 牛 亚 科 一 天 牛 亚 科 一 沟 腔 天 牛 亚 科 ,该 观点 
与 Linsley( 1961) 相似 ,不 同 的 是 异 天 牛 亚 科比 锯 天 
牛 亚 科 更 为 原始 ,着 天 牛 仍 作为 亚 科 存在 。Wang 和 
Chiang(1991 ) 结 合十 生物 学 与 十 地 理学 对 天 牛 各 亚 
科 的 起 源 问题 进行 了 人 研究 ,认为 沟 胜 天 牛 亚 科 最 为 
原始 ,其 分 化 顺 友 为 沟 腔 天 牛 亚 科 一 锯 天 牛 亚 科 一 
天 和 牛 亚 科 一 花 天 牛 亚 科 一 幽 天 牛 亚 科 或 者 为 沟 肥 天 
牛 亚 科 一 天 牛 亚 科 一 锯 天 牛 亚 科 一 花 天 牛 亚 科 一 网 
天 牛 亚 科 , 其 中 锯 天 牛 亚 科 包 括 异 天 牛 类 群 , 幽 天 牛 
亚 科 包括 椎 天 牛 类 群 , 花 天 和牛 亚 科 包括 膜 花 天 牛 类 
群 ; 王 文 遍 和 瑚 新 明 (2002 ) 以 天 牛 的 雌性 生殖 天 为 
寺 征 人 研究 表明 , 锯 天 牛 和 狭 胸 天 和 牛 的 分 类 地 位 不 确 
XE. ,演化 顺序 为 着 天 牛 科 一 椎 天 牛 亚 科 一 幽 天 牛 亚 
科 一 花 天 牛 亚 科 一 天 牛 亚 科 一 沟 肥 天 牛 亚 科 。 狭 胸 
天 牛 亚 科 的 归属 也 引起 广泛 的 争论 :Napp(1994 ) 并 
未 将 狭 胸 天 牛 列 人 天 牛 科 中 ,Svacha 等 (1997 ) 将 狭 
胸 天 牛 亚 科 归 为 暗 天 牛 科 ; 国 内 学 者 尹 新 明 (1999 ) 
的 研究 表明 狭 胸 天 牛 亚 科 和 上 暗 天 牛 科 的 大 部 分 特征 
较为 近 绿 ; 王 文 孔 和 尹 新 明 (2002 ) 的 人 研究 表明 锯 天 
牛 亚 科 和 狭 胸 天 牛 是 古老 类 群 ,但 是 其 归属 地 位 仍 
需要 进一步 研究 。 花 天 牛 亚 科 的 分 类 地 位 也 引起 广 
泛 的 讨论 :Crowson(1955 ) 认 为 花 天 牛 涉及 的 类 群 很 
T ,并 根据 形态 学 研究 认为 其 具有 古老 类 和 群 的 特征 ， 
然而 其 访 花 性 等 特征 应 该 和 地 球 上 的 显 花 植物 有 
关 , 因 此 其 演化 地 位 有 可 能 不 是 很 十 老 的 类 和 群 ,需要 
进一步 研究 。 
分 子 系统 学 的 人 研究 可 与 传统 的 形态 学 人 研究 相 
结合 ,能 为 解决 这 些 争 议 问题 提供 有 力 证 据 。 核 









































糖 体 18S rDNA 和 28S rDNA 由 于 其 保守 性 强 , 且 
携 刘 大 量 的 遗传 信息 ,已 经 在 鞘翅 目 昆 虫 的 系统 
进化 研究 上 得 到 了 应 用 (Hunt et al., 2007; 
Marvaldi et al., 2009) 。 然 而 关于 天 和 牛 高 阶 元 的 人 研 
REL, Hunt 等 (2007 ) 利用 核糖 体 18S rDNA 以 
及 线粒体 COL 和 16S rDNA 对 鞘翅 目的 高 级 阶 元 
系统 发 育 进行 了 人 研究 ,该 研究 中 涉及 到 了 天 牛 类 
群 ,其 中 异 天 牛 亚 科 、 饥 天 牛 亚 科 、 天 牛 亚 科 和 沟 
肥 天 牛 亚 科 位 于 一 文 , 作 为 天 牛 科 最 先 分 化 的 类 
群 , 异 天 牛 最 和 完 分 化 ,其 次 为 锯 天 牛 亚 科 、 沟 用 天 
牛 亚 科 与 天 牛 亚 科 成 姊妹 群 , 花 天 牛 亚 科 、 膜 花 天 
牛 亚 科 和 瘦 天 牛 亚 科 聚 为 一 文 , 暗 天 牛 亚 科 和 椎 
天 和 牛 亚 科 与 叶 甲 类 群 较 为 近 缘 。 我 国 利 用 分 子 系 
统 学 方法 对 天 牛 的 高 阶 元 研究 较 少 ( 张 健 等 ， 
2010,2012) 。 利 用 18S rDNA 与 28S rDNA 序列 联 
合 分 析 天 和 牛 高 级 阶 元 间 系 统 发 育 关 系 的 研究 尚未 
见报 道 。 本 研究 通过 18S rDNA V4 和 V7 区 以 及 
28S rDNA D2 和 D3 区 进行 测序 分 析 , 其 中 部 分 
18S rDNA 的 序列 来 源 于 魏 子 涵 等 (2014 ) 的 研究 
中 ,对 天 牛 总 科 2 BE 6 亚 科 共 63 种 ( 亚 种 ) 天 牛 进 
行 分 子 进化 研究 ,从 分 子 水 平 上 讨论 天 牛 亚 科 间 
的 进化 关系 ,并 探讨 核糖 体 DNA. 联合 分 析 方 法 是 
否 适合 于 天 牛 高 阶 元 的 人 研究 ,为 天 牛 总 科 的 系统 
发 育 人 研究 提供 基础 依据 。 
































1 材料 与 方法 


1.1 供 试 虫 源 

供 试 虫 样 及 详细 来 源 见 表 1 ,研究 涉及 的 63 种 
天 牛 均 为 成 虫 。 标 本 采集 时 间 为 2008 -2012 年 , 主 
要 集中 在 每 年 的 6 -10 月份。 采集 地 点 主要 为 河南 
省 ,少数 采集 于 广州 南 岭 和 贵州 等 地 。 和 采集 到 的 标 
本 用 无 水 乙醇 浸泡 , 其 中 单 锥 育 天 和 牛 Thranius 
simplex Gahan 和 双 市 粒 示 天 牛 Lamiomimus gottschei 
Kolbe 为 干 标 本 。 外 群 选 择 为 河南 农业 大 学 标本 馆 
提供 的 肖 叶 甲 科 ( Eumolpidae) f9 4E IEEE i np IE 
Chrysochus chinensis Baly 以 及 GenBank 数据 库 中 象 
甲 科 ( Curculionidae ) 的 Araucarius minor Kuschel 
( GenBank 序列 号 : FJ867736 和 AF308351), 

KWARK ^F pi A AE MEUBLE Lamiinae ) 
28 种 ,天 牛 亚 科 (Cerambycinae ) 24 种 , 锯 天 牛 亚 科 
( Prioninae) 3 种 , 花 天 牛 亚 科 (Lepturinae) 2 f, Ki 
天 牛 亚 科 (Aseminae) 3 种 , 狭 胸 天 牛 亚 科 (Philinae ) 
2 种 , 瘦 天 牛 科 (Disteniidae) 1 种 ( 表 1)。 
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R1 标本 采集 信息 及 18S rDNA V4 $0 V7 区 以 及 28S rDNA D2 和 D3 区 序列 登录 号 
Table 1 Specimens collecting details and GenBank accession numbers of 18S rDNA ( V4 and V7 region) 
and 288 rDNA (D2 and D3 region) sequences 
GenBank 登录 号 
FEE : 采集 地 采集 时 间 GenBank accession no. 
: ; 种 类 Species . . 
Family/ subfamily Collecting locality Collecting date 18S rDNA 28S rDNA 
V4 V7 D2D3 
BERF 河南 平顶山 
2012.6 KF141946 KF142010 KF142075 
Psacothea hilaris (Pascoe) Pingdingshan, Henan 
JE/H ECKE 河南 平顶山 
2012.0 KF141948 KF142012 KF142077 
Anoplophora glabripennis ( Motschulsky ) Pingdingshan, Henan 
黄斑 星 天 牛 河南 平顶山 
jn 2008. 6 KF141947  KF142011 KF142076 
Anoplophora nobilis ( Ganglbauer ) Pingdingshan, Henan 
星 天 牛 河南 台州 
2010.7 KF141949 KF142013  KF14207 
Anoplophora chinensis ( Forster) Yuzhou, Henan s 
ERF 河南 平顶山 
2011.6 KF141950  KF142014  KF142079 
Anoplophora lurida ( Pascoe) Pingdingshan, Henan 
ary y XT 
FRERE l 河南 平顶山 2011.6 KF141951 — KF142015  KF142080 
Uraecha angusta. ( Pascoe) Pingdingshan, Henan 
HE y zn| y AM 
: BUPRGSUAAE 河南 贺州 2010 KF141953 KF142017  KF142082 
Lamiomimus gottschei Kolbe Yuzhou, Henan 
ENERE 河南 平顶山 
2011.6 KF141952  KF142016 KF142081 
Paraglenea fortunei ( Saunders) Pingdingshan, Henan 
A DAK N 住 
FARER : 河南 焦作 2012.5 KF141954  KF142018  KF142083 
Glenea sp. Jiaozuo, Henan 
暗 翅 简 天 牛 河南 南阳 
2010 KF141955  KF142019 KF142084 
Oberea fuscipennis ( Chevrolat ) Nanyang, Henan 
问 大 咎 属 . 河南 焦作 2012.5 KF141956 KF142020 KF142085 
Oberea sp. 1 Jiaozuo, Henan 
简 河南 焦 
问 大 牛 属 河南 焦作 2012.5 KF141957 KF142021  KF142086 
Oberea sp. 2 Jiaozuo, Henan 
ile : 三 条 小 简 河南 焦 
EREE 三 系 小 简 大 和 牛 l 河南 焦作 2012.5 KF141958 KF142022 KF142087 
Lamiinae Phytoecia sibirica ( Gebler) Jiaozuo, Henan 
双 簇 污 天 牛 河南 济源 
2012.5 KF141959  KF142003  KF142088 
Moechotypa diphysis ( Pascoe) Jiaozuo, Henan 
—BREH ERE 广东 南 岭 
2012.5 KF141960 KF142004  KF1420; 
Xenohammus bimaculatus Schwarzer Nanling, Guangdong d 
贵州 茂 兰 自然 保护 
PITE 贯 州 友 兰 日 然 你 护 区 
Maolan National Natural 2012.5 KF141961 — KF142025  KF142090 
Monochamus alternatus Hope . 
Reserve, Guizhou 
E Ju ETE 
RAAR E . : 河南 平顶山 2011.7 KF141962  KF142026  KF142091 
Ropica subnotata Pic Pingdingshan, Henan 
柳 坡 天 牛 河南 平顶山 
2012.6 KF141963 KF142027  KF142092 
Pterolophia rigida ( Bates ) Pinegdingshan, Henan e: 
贵州 成 兰 自然 保护 区 
BISHER: bis a i 
: : Maolan National Natural 2012.5 KF141964  KF142028  KF142093 
Phelipara radowskyi Hua 
Reserve, Guizhou 
M 贵州 茂 兰 自然 保护 区 
RARE 
f f Maolan National Natural 2012.5 KF141965 KF142029 KF142094 
Exocentrus guttulatus subconjunctus ( Gressitt ) 
Reserve, Guizhou 
HEE Sz E ARAE 
Maolan National Natural 2012.5 KF141966 KF142031 KF142096 
Olenecamptus clarus Pascoe : 
Reserve, Guizhou 
中 华 八 星 粉 天 牛 河南 南阳 
2009.7 KF141968 KF142032  KF142 
Olenecamptus octopustulatus chinensis Dillion & Dillion Nanyang, Henan i 
一 一 了 A. y NA In 
云 更 日 条 大 午 河南 平顶山 2012.0 KF141969 KF142033  KF142098 


Batocera lineolata Chevrolat 


Pingdingshan, Henan 
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续 表 1 Table 1 continued 
GenBank 登录 号 
科 亚 科 采集 地 采集 时 间 GenBank accession no. 
i l 种 类 Species l 
Family/ subfamily Collecting locality Collecting date 18S rDNA 28S rDNA 
V4 V7 D2D3 
xe y XT 
RRRA ; 河南 平顶山 2012.6 KF141970  KF142034  KF142099 
Apriona germari (Hope) Pingdingshan, Henan 
ERE EAR > 
; . Maolan National Natural 2012.5 KF141972 | KF142036  KF142101 
Mesosa atrostigma. Gressitt ] 
Reserve, Guizhou 
占 y à 小 上 
VASE 河南 潮州 2011.8 KF141973 KF142037 KF142102 
Mesosa myops ( Dalman) Zhengzhou, Henan 
AN y XT 
BARKE À . 河南 平顶山 2012.6 KF141974 KF142038  KF142103 
Mesosa irrorata Gressitt Pingdingshan, Henan 
WERE 河南 郑州 
2012.6 KFI41971  KF142035  KF1421 
Pseudanaesthetis langana Pic Zhengzhou, Henan ui 
B X < IA 
HIRE 三 东南 * 2012. 5 KF141967 KF142030 KF142095 
Trachylophus sinensis Gahan Nanling, Guangdong 
Zi X YTA 
RERE . 河南 平顶山 2011.6 KF141975  KF142039 KF142104 
Xylotrechus pyrrhoderus Bates Pingdingshan, Henan 
EA y XT 
TERCERA 河南 平顶山 2012.6 KF141976 KF142040 KF142105 
Rhaphuma elongata Gressitt Pingdingshan, Henan 
六 斑 绿 虎 天 牛 河南 平顶山 
2012.6 KF141977 — KF142041  KF142106 
Chlorophorus sexmaculatus ( Motschulsky ) Pingdingshan, Henan 
河南 焦 
KEJE ; 河南 焦作 2012.5 KF141978 KF142042 KF142107 
Clytus sp. 2 Jiaozuo, Henan 
EAK X fi 
ESI IR DEAE LEE 河南 焦作 2012.5 KF141979  KF142043  KF142108 
Teratoclytus plavilstshikovi Zaitz Jiaozuo, Henan 
弧 纹 绿 虎 天 牛 河南 郑州 
2012.6 KF141980  KF142044  KF1421 
Chlorophorus miwai Gressitt Zhengzhou, Henan Ne 
cp y XT 
ERKE . 河南 平顶山 2011.6 KF141981 KF142045  KF142110 
Clytus sp.1 Pingdingshan, Henan 
KERF 河南 郑州 
2011.6 KF141982 — KF142046  KF142111 
Trichoferus campestris ( Faldermann) Zhengzhou, Henan 
HH Jürg VT 
RURE ; l 河南 平顶山 2012.6 KF141983  KF142047  KF142112 
xS Massicus raddei ( Blessig) Pingdingshan, Henan 
MES 
. | "H ME 页 
Cerambycinae RARE ) PFD Lu 2012.5 KF141984  KF142048  KF142113 
Trirachys orientalis Hope Pingdingshan, Henan 
SL 河南 郑州 
REAR IRIAN 2011.6 KF141985 KF142049 KF142114 
Margites fulvidus (Pascoe) Zhengzhou, Henan 
LEY 2 y SZ 页 
"a A 家 天 咎 | . 河南 下 本 "i 2011.6 KF141986 — KF142050  KF142115 
Zegriades aurovirgatus Gressitt Pingdingshan, Henan 
HESS 河南 焦 
BEKAK mt 河南 焦作 2012.5 KF141987 KF142051  KF142116 
Purpuricenus sideriger Fairmaire Jiaozuo, Henan 
IEEE 贵州 茂 兰 自然 保护 区 
dado Maolan National Natural 2012.5 KF141988 — KF142052  KF142117 
Purpuricenus temminckii ( Guerin-Meneville) 
Reserve, Guizhou 
东 2n Jn YE 页 
东亚 SAVE 河南 下 本 加 2012.6 KF141989  KF142053  KF142118 
Allotraeus asiaticus (Schwarzer) Pingdingshan, Henan 
TARE 河南 平顶山 
2012.6 KF141990  KF142054  KFI4211 
Nadezhdiella cantori ( Hope) Pingdingshan, Henan t 
i 河南 焦 
AEGRIS AE — 河南 焦作 2012.5 KF141991 KF142055  KF142120 
Pyrestes rugicollis Fairmaire Jiaozuo, Henan 
A o A. Y H M^. In 
全 欢 双 条 入 千 n 河南 平顶山 2012.6 KF141992 KF142056  KF142121 
Xystrocera globosa (Olivier) Pingdingshan, Henan 
` RIN "H XT 
UNA 河南 平顶山 2012.6 KF141993 KF142057 KF142122 


Stenygrinum quadrinotatum Bates 


Pingdingshan, Henan 
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续 表 1 Table 1 continued 
GenBank 登录 号 
科 亚 科 采集 地 采集 时 间 GenBank accession no. 
i l 种 类 Species l 
Family/ subfamily Collecting locality Collecting date 18S rDNA 28S rDNA 
V4 V7 D2D3 
FRE 河南 平顶山 
2012.0 KF141994 KF142058 KF142123 
Rhytidodera bowringii White Pingdingshan, Henan 
ydi N MZ 而 
"her AE . 河南 平顶山 2009.7 KF1419905  KF142059 KF142124 
Thranius simplex Gahan Pingdingshan, Henan 
HIIRE 河南 郑州 
2011.6 KF141996  KF142060 KF142125 
Aromia bungii ( Faldermann) Zhengzhou, Henan 
Azt N YE 
Zu Xt 河南 平顶山 
2011.7 KF141997 KF142061  KF142126 
Polyzonus fasciatus (Fabricius) Pingdingshan, Henan 
亦 短 梗 天 牛 广东 广州 
2011.7 KF141998 — KF142062 KF142127 
Arhopalus unicolor ( Gahan) Guangzhou, Guangdong 
A N 1 y Uhh 
BORIEERI ERE . 河南 郑州 2011.7 KF141999 KF142063 KF142128 
Aseminae Arhopalus rusticus ( Linnaeus) Zhengzhou, Henan 
CJ 新 疆 秦 岭 
i XCTI A. . HER ;! 2012. 6 KF142002  KF142066  KF142131 
Tetropium castaneum ( Linnaeus) Qinling, Xinjiang 
: FERA E 三 东 珠 海 2012.4 KF142001  KF142065  KF142130 
狭 胸 天 和 牛 亚 科 Philus antennatus ( Gyllenhal) Zhuhai , Guangdong 
Philinae TER RIR 
. MIMIKE 三 东 珠 海 2012.4 KF142000 KF142064 KF142129 
Philus pallescens Bates Zhuhai, Guangdong 
中 华 薄 翅 天 牛 河南 郑州 
2012.6 KF142003 KF142067 KF142132 
Megopis sinica ( White) Zhengzhou, Henan 
ENA melo TIRAR 2010 KF142004 KF142068 KF142133 
Prioninae Dorysthenes fossatus Pascoe Yuzhou, Henan 
; HARRIKA 三 东 珠 海 2012.5 KF142005  KF142069  KF142134 
Priotyrranus closteroides (Thomson ) Zhuhai, Guangdong 
Je Bj 河南 焦 
瘙 胸 金 花 大 千 : 河南 焦作 2012.5 KF142006  KF142070  KF142135 
花 天 牛 亚 科 Gaurotes tuberculicollis ( Blanchard) Jiaozuo, Henan 
Lepturinae } Zk pal 
: PIDIERA : 三 东南 2012.5 KF142007  KF142071  KF142136 
Leptura meridiosinica Gressitt Nanling, Guangdong 
N Zi N ATP 
EKER IURI. . : 河南 平顶山 2011.7 KF142008 KF14207/2  KF142137 
Disteniidae Cyrtonops asahinai Mitono Pingdingshan, Henan 
外 群 Outgroup 
肖 叶 甲 科 HEATERS Pap HH 河南 新 密 
2009 KF142009 KF142073  KF142074 
Eumolpidae Chrysochus chinensis Baly Xinmi, Henan 
ARE Araucarius minor Kuschel FJ867736 | FJ867736 AF308351 


Curculionidae 


1.2 基因 组 DNA 的 提取 

取 天 牛 胸 部 肌肉 约 25 mg 于 1.5 mL 的 离心 管 
里 ,加 入 250 pL ARRIA 20 pL EAR K 上海 
生 工 ) ,人 研磨 至 天 牛 的 胸部 肌肉 组 织 完全 裂解 ,随后 步 
又 根据 基因 组 试剂 盒 ( 上 海 莱 枫 ) 说明书 进行 操作 。 
利用 NanoDrop 1000 X 5h 4j 36 26 JÆ il ( Thermo 
Scientific, America) 对 所 抽 提 出 的 总 DNA 进行 浓度 
IAE BE ( OD4,/0D», ) 检测 。 将 检测 大 于 100 ng/ pL 
基因 组 DNA 稀释 至 80 ng/L 左右 ,并 分 为 两 份 ,一 
份 用 于 PCR 扩 增 ,一 份 保存 于 -20%C 冰 箱 。 
1.3 目的 片段 扩 增 

利用 4 对 引物 ( 表 2) 分 别 对 18S rDNA V4 和 





V7 区 以 及 28S rDNA D2 和 D3 区 序列 进行 扩 增 ,其 
中 18S rDNA V4 和 V7 区 都 在 300 ~ 400 bp, 28S 
rDNA D2 区 为 500 bp 左右 ,28S rDNA D3 区 为 200 ~ 
300 bp, PCR 反应 采用 25 uL 体系 ,包括 12.5 uL 
的 2 x Taq PCR Master Mix ( E, JEM) , 正 反 引物 
各 1 uL,9. 5 uL ddH,0,1 pL 的 模板 DNA, 28S 
rDNA D2 和 D3 区 的 扩 增 条 件 为 :94% 预 变 性 3 min 
后 , 按 以 下 参数 进行 33 个 循环 ,94Y 变性 40 s,55*C 
退火 1 min,72% 延伸 1 min, fff Y SE tE Je; 72°C 延伸 
10 min, 18S rDNA V4 和 V7 区 的 PCR 扩 增 ,退火 温 
度 为 67%C ,其 他 步骤 与 28S rDNA D2 和 D3 区 相同 。 
扩 增 产物 用 1.5% BN Bs REPRE E DK E ÍT Tr 0 , 对 
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扩 增 结果 展 好 的 样品 通过 PCR 纯化 试剂 盒 纯 化 
( Axygen Scientific ,杭州 ) ,之 后 委托 上 海 生 工 公司 





对 样品 进行 正 反 双向 测序 , 测序 仪 为 ABL 
PRISM3730 ,测序 试剂 为 BigDyeterminator v3. 1, 


表 2 本 研究 所 用 的 引物 
Table2 Primers used in this study 








引物 名 称 引物 序列 (5 -3) 目的 片段 参考 文献 
Primer name Primer sequence Target fragment References 
V4F TGGTGCCAGCAGCCGCGGTA A 
188 rDNA V4 Raupach et al., 2010 
V4R CCTCTAACGTCGCAATACGAATGCCC 
V7F CTTAAAGGAATTGACGGAAGGGCACCACC 
18S rDNA V7 Raupach et al., 2010 
V7R GATTCCTICAGTGTAGCGCGCGTG 
D2F GAGTTCAAGAGTACGTGAAACCG 
288 rDNA D2 Gillespie et al., 2003 
D2R CCTTGGTCCGTGTTTCAAGAC 
D3F GGACCCGTCTTGA AACAC 
28S rDNA D3 Raupach et al., 2010 
D3R GCATAGTTCACCATCTTTC 


1.4 DNA 目的 序列 处 理 及 同 质 性 测验 (incongruence 
length difference test, ILD) 

将 获得 的 序列 分 别提 交 至 NCBI( http ;//blast. 
ncbi. nlm. nih. gov/Blast. egi) 进行 同 源 性 比 对 以 确 
保 序 列 的 准确 性 。 用 Chromas 2 (Technelysium Pty 
Ltd., Australia) 对 测序 获得 的 DNA 序列 进行 峰值 观 
邃 ,并 用 软件 MEGA 5. 1 ( Tamura et al., 2011) 进行 
正 反 链 拼接 ,最 后 将 每 个 样本 4 个 区 的 DNA 序列 进 
行 拼接 并 比 对 。 由 于 28S rDNA 和 18S rDNA 具有 
不 同 的 遗传 速率 ,因此 在 进行 联合 序列 分 析 前 利用 
PAUP -4. Ob10 ( Swofford, 2002) 进行 同 质 性 检验 
( Farris et al., 1994) ,来 判断 数据 是 否 可 以 联合 
分 析 。 

1.5 序列 组 成 分 析 及 遗传 距离 分 析 

经 过 比 对 后 ,拼接 好 的 序列 全 长 为 1 404 bp( 包 
含 gap) ,用 MEGA 5.1 统计 核 音 酸 A,T,C,G 的 组 成 、 
保守 位 点 (conserved sites) 比例 、 转 换 值 (transition , 
Ts) HRR (transversion, Tv) ,并 计算 亚 科 间 的 遗传 
距离 ,距离 模型 选用 非 矫 正 p 距离 (uncorrected p- 
distance) ,并 对 基因 饱和 上 度 进行 分 析 。 

1.6 系统 发 育 树 的 构建 

分 别 用 邻接 法 (NJ) .最 大 似 然 法 (ML) 和 贝 叶 
斯 推论 法 ( BI) 构建 进化 树 , 其 中 ML 法 和 BI 法 利用 
Modeltest 3. 7 ( Posada and Crandall, 1998) 计算 进化 
模型 。 联 合 序 列 最 佳 蔡 换 模 型 为 general time 


reversible( GTR) + the proportion of invariable sites 





(I) + the gamma-shaped distribution of rates across 
sites( G) 。 最 大 似 然 法 (ML) 建树 由 在 线 生 物 信息 
学 工具 phyML 3. 0 ( Guindon et al., 2005) 构建 
( http :// www. atgc-montpellier. fr/phyml/) ,参数 设置 
如 下 :替代 模型 (substitution model ) 选择 “GTR”, 不 


变 位 点 比例 (proportion of invariable sites ) 选择 
“estimated”, 替代 率 类 别 数 (number of substitution 
rate categories) 设 为 "6”, 树 形 提高 策略 选择 “SPR ” ， 
同时 运行 SH-like 测试 对 进化 树 评 佑 ,其 余 参数 用 默 
认 设 置 。 建 树 的 结果 用 FigTree 1.3.1(Rambaut and 
Drummond, 2009) 软件 打开 。 

运用 MrBayes3. 2. 1 ( Ronquist and Huelsenbeck , 
2003) 构建 中叶 斯 树 , 同 时 运行 两 个 独立 的 分 析 ， 
个 分 析 包 含 4 条 马尔 科 夫 链 (3 条 热 链 ,1 条 冷 链 ) ， 
分 析 从 随机 树 起 始 , 共 运行 5 000 000 代 , 每 500 代 
取 一 次 样 ,将 前 25% 样 本 作为 老化 样本 舍 去 ,余下 
的 用 作 构 建 进化 树 和 统计 后 验 概率 (posterior 
probability, PP) , 当 5 000 000 代 运 行 完成 ,分 裂 频 
率 标准 差 ( average standard deviation of the split 
frequency) RF 0. 01 时 ,认为 两 个 独立 分 析 的 结 
非常 接近 。 由 于 该 分 析 是 联合 基因 序列 分 析 , 因 此 
采用 划分 区 段 方式 建树 ,分 别 对 两 个 基因 序列 进行 
模型 检测 ,通过 Modeltest3. 7 ( Posada and Crandall, 
1998 ) 软件 分 析 得 出 ,2 个 基因 序列 的 最 佳 模 型 均 为 
GTR + I + CC, 建树 结果 用 FigTree 1. 3. 1 ( Rambaut 
and Drummond 2009) 进行 查看 。 








2 结果 与 分 析 


2.1 28S rDNA D2 和 D3 区 以 及 18S rDNA V4 和 
V7 区 序列 的 测定 与 基因 登录 号 

对 28S rDNA D2 和 D3 区 以 及 18S rDNA V4 和 
V7 区 序列 测序 ,涉及 到 63 种 天 牛 。 利 用 4 对 引物 
分 别 扩 增 ,结果 显示 ,18S rDNA 的 V4 和 V7 区 都 在 
300 ~ 400 bp 左右 ,28S rDNA D2 区 为 500 bp 左右 ， 
D3 区 为 200 ~300 bp, 符 合理 论 结果 。63 种 天 牛 基 
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因 登 录 号 18S rDNA V4, V7: KF141946 ~ KF142072 ; 
28S rDNA D2D3 :KF142075 ~ KF142137 (2 1) 。 
2.2. 基因 的 同 质 性 检验 及 保守 位 点 分 析 

利用 PAUP 4. ObIO ( Swofford, 2002) 对 序列 进 
行 同 质 性 检验 ,结果 28S rDNA 和 18S rDNA 序列 可 
以 进行 联合 分 析 (P = 0. 103 ) 。 基 因 拼 接 并 进行 同 
源 性 比 对 后 序列 总 长 为 1 404 bp。 运 用 MEGA 5.1 
对 核 童 酸 组 成 的 分 析 , 结 采 表明 : 碱 基 A, T, C LG 
的 平均 含量 分 别 为 19. 796, 23. 3% ，26. 2% 和 
30.8% ,可 变 位 点 446 个 , 占 总 位 点 的 32. 0% ,保守 
位 点 958 个 , 占 总 位 点 的 68. 0% ,转换 / 题 换 平 均值 








(R 值 ) 为 1.73 ,转换 值 高 于 题 换 值 。 
2.3 遗传 距离 分 析 

利用 MEGA 5.1 软件 对 各 个 亚 科 间 进行 非 矫正 
p 距离 的 计算 ( 表 3) , 结 采 显示 :外 群 与 天 牛 类 和 群 各 
亚 科 的 遗传 距离 均 大 于 天 牛 类 群 各 亚 科 之 间 的 遗传 
距离 , 饮 天 牛 亚 科 与 从 天 牛 亚 科 距离 最 小 (0. 040 ) 。 
与 沟 肥 天 牛 亚 科 距离 最 近 的 是 天 牛 亚 科 (0. 059), 
其 次 为 花 天 牛 亚 科 (0. 061) ;与 天 牛 亚 科 距离 最 近 
的 是 饮 天 牛 亚 科 (0. 043 ) , 其 次 为 花 天 牛 亚 科 
(0.045) 。 与 幽 天 牛 亚 科 距离 最 近 的 为 天 牛 亚 科 和 
花 天 牛 亚 科 (0.063) 


RI 亚 科 之 间 非 矫正 p 距离 (下 三 角 ) 以 及 标准 差 ( 上 三 角 ) 


Table 3 Uncorrected p-distance (lower triangle) and standard deviation (upper triangle) between the subfamilies 





PEE Family/subfamily N 1 3 4 5 6 7 8 

1 VIE AE ERE Laminae 28 0.004 0.005 0.005 0.005 0.005 0.006 € 0.006 
2 天 牛 亚 科 Cerambycinae 24 0.059 0.004 0.004 0.005 0.005 0.006 € 0.006 
3 花 天 牛 亚 科 Lepturinae 2 0.061 0.045 0.004 0.006 0.005 0.007 0.006 
4 锯 天 牛 亚 科 Prioninae 3 0.062 0.043 0.040 0.006 0.005 0.007 0. 006 
5 幽 天 牛 亚 科 Aseminae 3 0.072 0.063 0.063 0.066 0.006 0.007 0.006 
6 狭 胸 天 和 牛 亚 科 Philinae 2 0.071 | 0.056 | 0.055 0.059 . 0.068 0.007 . 0.006 
7 瘦 天 牛 科 Disteniidae 1 0.081 0.077 0.077 0.079 0.076 0.083 0.007 
8 外 群 Outgroup 2 0. 093 1.000 0.084 0.084 0.098 0.093 0.100 


N: 物种 数 Number of species. 


2.4 碱 基 置换 饱和 度 分 析 

以 非 校正 p 距离 为 村 坐标, 以 各 个 碱 基 的 转换 
(Ts) 和 匡 换 (Tv) 为 纵 坐 标 作 散 点 图 ,进行 碱 基 蔡 换 
饱和 度 分 析 。28S rDNA 和 18S rDNA 以 及 联合 序列 
的 饱和 度 分 析 绪 有 末 ( 图 1 ) 显示 : 碱 基 置 换 数 与 非 校 
正二 距离 呈 线 性 关系 , 碱 基 突变 未 达到 饱和 ,可 以 用 
于 系统 进化 树 分 析 。 
2.5 ”进化 树 拓扑 结构 分 析 

以 NJ 法 .ML 法 和 BI 法 3 种 方法 建树 ,结果 见 图 
2 ~4。 不 同方 法 构建 进化 树 均 得 出 相似 的 拓扑 结构 
上 且 具 有 较 高 的 置信 和 度 。 进 化 树 的 结果 表明 : 供 试 天 牛 
类 群 形成 一 文 , 与 外 群 明显 分 开 , 且 具有 较 高 的 置信 
度 :NJ 树 bootstrap 值 为 81( 图 2),ML 树 的 SH-Like 
值 为 84( 图 3) ,BI 树 的 的 后 验 概率 (PP) 值 为 89。 建 树 
结果 均 显 示 天 牛 类 群 主 要 分 化 出 3 个 进化 文 (A,B 和 
C) :最 完 分 化 出 的 A 校 由 瘦 天 牛 科 和 幽 天 牛 亚 科 组 成 ; 
然后 进化 树 分 化 出 了 B AA C 校 ,其 中 B 枝 由 锯 天 牛 
亚 科 ` 狭 胸 天 牛 亚 科 、 花 天 牛 亚 科 和 天 牛 亚 科 组 成 ,C 
KIARA MEE, 3 种 建树 方法 的 结果 都 支持 沟 肥 














天 牛 亚 科 天 牛 亚 科 和 人 花 天 牛 亚 科 为 单 系 进 化 群 。 





60 
45 
30 
15 


o aita 
0 0.05 





Number of substitutions 


0.10 0.15 0.20 


0.25 0.30 


120 : 
100 Combined se 





0 0.02 0.04 0.06 0.08 0.10 0.12 0.14 
pE p distance 
图 1 18S rDNA 28S rDNA 及 联合 序列 饱和 度 分 析 


Fig. 1 
28S rDNA and the combined sequence 
Ts: 转换 值 Transition; Tv; 颠 换 值 Transversion. 


Saturation degree analysis of 18S rDNA, 
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97 58 Oberea fuscipennis 
berea sp.1 
Oberea sp.2 
Phytoecia sibirica 
Paraglenea fortunei 
Glenea sp. 
63 Moechotypa diphysis 
Pseudanaesthetis langana 
— Exocentrus guttulatus subconjunctus 
~ 92 -Anoplophora nobilis 
100 56 Anoplophora glabripennis 
78 Anoplophora lurida 
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图 2 基于 18S rDNA VA 和 V7 区 以 及 28S rDNA D2 和 D3 区 构建 的 距离 树 进化 树 
Fig.2 NJ phylogeny tree based on 18S rDNA (V4 and V7 region) and 28S rDNA (D2 and D3 region) sequences 


A 枝 由 瘦 天 牛 科 和 幽 天 牛 亚 科 组 成 ; B 枝 由 狭 胸 天 牛 亚 科 、 锯 天 牛 亚 科 、 花 天 牛 亚 科 和 天 牛 亚 科 组 成 ; C TS ER AREER , o FN S 
位 核 苷 酸 的 蔚 换 量 ; 图 3 和 4 同 。 进 化 枝 显 示 大 于 50 的 bootstrap 值 。A branch includes Disteniidae and Aseminae; B branch includes Philinae, 


Prioninae ，Lepturinae and Cerambycinae; C branch includes Lamiinae. The scale represents expected number of substitutions per site. The same for Figs. 


3 and 4. Bootstrap value above 50 was labeled on the interior branch. 


3 讨论 


本 研究 成 功 扩 增 出 63 种 天 牛 18S rDNA V4 和 
V7 区 以 及 28S rDNA D2 和 D3 区 序列 。 经 MEGA 
5.1 软件 分 析 后 :18S rDNA V4 和 V7 区 序列 693 
bp ,保守 位 点 608 个 (88.0% ) ,28S rDNA D2 和 D3 
区 序列 711 bp, 保守 位 点 350 个 (49. 096) 。18S 
rDNA VA 和 V7 区 序列 比 28S rDNA D2 和 D3 区 序 
列 更 加 保守 ,这 与 刘 殿 锋 和 蒋 国 芳 (2005 ) 报道 的 一 
致 。 序 列 的 同 质 性 检验 结果 忆 > 0.05( Farris et al., 
1994) ,因此 可 以 进行 联合 。 将 18S rDNA 与 28S 
rDNA 序列 进行 联合 ,在 序列 饱和 度 分 析 中 ,联合 序 


列 分 析 比 单 基 因 序 列 分 析 具 有 更 好 的 线性 关系 (图 
1) ,因此 更 适合 于 进化 树 的 构建 ,利用 联合 序列 构 
建 的 进化 树 均 具有 相似 的 拓扑 结构 和 较 高 的 解 
析 度 。 

分 别 利用 NJ 法 、ML 法 和 BI 法 对 联合 序列 构 
建 进化 树 ,3 种 方法 建树 具有 相似 的 拓扑 结构 : 进 
化 树 基部 (A 枝 ) 由 瘦 天 牛 科 与 幽 天 牛 亚 科 组 成 暗 
示 其 较为 近 缘 ,然而 3 种 进化 树 的 置信 度 存 在 差 
异 较 大 ,其 中 NJ 树 中 的 bootstrap 值 为 67( 图 2)， 
ML 树 中 的 SH-Like 值 为 58( 图 3), 而 BI 树 中 的 
PP 值 为 92( 图 4), 男 外 未 校正 p 距离 显示 瘦 天 牛 
科 与 天 牛 亚 科 之 间 的 距离 (0.076 ~0.100) 大 于 天 
牛 各 亚 科 之 间 的 距离 (0.040 ~0.072) ,因此 瘦 天 
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图 3 基于 18S rDNAV4 fil V7 区 以 及 28S rDNA D2 和 D3 区 构建 的 最 大 似 然 进化 树 
Fig.3 ML phylogeny tree based on 18S rDNA (V4 and V7 region) and 28S rDNA (D2 and D3 region) sequences 
进化 枝 显 示 大 于 50 的 SH-Like 值 。SH-Like value above 50 was labeled on the interior branch. 


牛 是 否 独 立成 科 需 要 进一步 的 研究 。 进 化 树 中 部 
(B 枝 ) 由 天 牛 亚 科 、 饮 天 牛 亚 科 、 花 天 牛 亚 科 和 狭 
胸 天 牛 亚 科 组 成 。 狭 胸 天 牛 类 的 归属 一 直 存 在 争 
议 : 国 外 学 者 Švácha 等 (1997 ) 将 狭 胸 天 牛 作为 暗 
天 牛 科 ( Vesperidae ) 的 一 个 亚 科 , 而 我 国学 者 王 直 
诚 (2003 ) 和 华 立 中 等 (2009 ) 分 别 将 狭 胸 天 牛 归 
为 饥 天 牛 亚 科 的 一 个 属 。 本 研究 结果 表明 : 狭 胸 
天 牛 位 于 B 枝 的 基部 ,最 早 分 化 出 来 ,然后 分 化 为 
锯 天 牛 亚 科 和 花 天 牛 亚 科 。 狭 胸 天 牛 应 该 为 较为 
古老 的 类 群 ,这 与 尹 新 明 (1999 ) 的 观点 相同 ,但 是 
狭 胸 天 牛 并 没有 出 现在 锯 天 牛 亚 科 的 分 支 中 , 因 
此 狭 胸 天 牛 应 当 独 立 作 为 一 个 亚 科 从 锯 天 牛 亚 科 
中 分 离 出 来 并 暂时 归 为 天 牛 科 。B 枝 中 句 天 牛 亚 
科 的 单 系 性 并 没有 得 到 一 致 的 结论 ,其 中 ML 树 
(图 3) 与 BI 树 (图 4) 均 文 持 了 饮 天 牛 亚 科 的 单 系 
性 ,而 NJ 树 (图 2) 并 没有 支持 该 结论 , 锯 天 牛 亚 科 
中 的 中 华 薄 姻 天 牛 M. sinica 最 先 分 化 但 其 
bootstrap 值 (56 ) 较 低 。 除 此 之 外 ,进化 树 都 表明 

















饮 天 牛 亚 科 和 花 天 牛 亚 科 先 聚 在 一 起 ,随后 又 与 
天 牛 亚 科 聚 在 一 起 ,该 结果 与 张 健 等 (2012 ) 利用 
28S rDNA 的 研究 相似 。 花 天 牛 的 分 类 地 位 在 本 人 研 
究 存 在 不 确定 性 , 王 齐 (1986) 以 十 植物 学 和 地 理 
学 进行 研究 认为 花 天 牛 的 访 花 性 以 及 大 多 有 昼 出 
性 ,这 与 显 花 植物 的 出 现 有 关 , 因 此 不 应 是 古老 的 
类 和 群 。 将 书 楠 和 陈 力 (2001 ) 对 花 天 牛 的 地 位 进行 
了 探讨 ,认为 其 应 该 处 于 幽 天 牛 亚 科 和 天 牛 亚 科 
之 间 较 进化 的 类 群 ,本 研究 中 花 天 牛 处 于 幽 天 牛 
亚 科 与 天 牛 亚 科 之 间 , 文 择 了 该 结论 。C 术 由 沟 肥 
天 牛 亚 科 构成 为 单 系 群 ,与 了 B CCS BEDS AS ,分 
化 晚 于 首 天 牛 科 和 了 幽 天 牛 亚 科 , 虽 然 其 进化 性 还 
需 进 一 步 人 研究 ,然而 本 研究 结果 已 与 Wang 和 
Chiang(1991 ) 等 观点 不 同 。 

本 实验 利用 18S rDNA V4 和 V7 区 以 及 28S 
rDNA D2 和 D3 区 的 研究 表明 :联合 分 析 在 天 牛 类 
群 亚 科 内 以 及 亚 科 间 都 有 较 高 的 解析 度 ,可 以 作为 
天 牛 高 级 阶 元 研究 有 效 的 分 子 标记 ;在 天 牛 高 级 阶 
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图 4 分 别 对 18S rDNA V4 和 V7 区 以 及 28S rDNA D2 和 D3 区 分 区 段 基因 定义 模型 构建 的 贝 叶 斯 树 
Fig.4 BI phylogeny tree based on 18S rDNA (V4 and V7 region) and 28S rDNA (D2 and D3 region) sequences 
by using partitioned scheme posterior probability 
进化 枝 显 示 大 于 50 的 后 验 概率 值 。PP value above 50 was labeled on the interior branch. 


元 的 研究 中 我 们 提供 了 一 些 有 用 的 信息 ,比如 证 明 
了 沟 肥 天 牛 亚 科 、 天 牛 亚 科 和 花 天 牛 亚 科 的 单 系 性 ， 
狭 胸 天 牛 的 独立 成 为 亚 科 等 。 然 而 目前 核 基因 联合 
分 析 的 局 限 是 ,在 天 牛 类 群 中 核 基因 18S rDNA V4 
和 V7 区 以 及 28S rDNA D2 和 D3 区 序列 在 GenBank 
数据 库 中 同 源 序列 并 不 是 很 多 ,因此 对 天 牛 总 科 高 
级 阶 元 进化 研究 会 产生 影响 ,比如 本 实验 中 缺少 上 暗 
天 牛 的 序列 ,无 法 去 验证 狭 胸 天 牛 亚 科 是 否 应 该 归 
为 暗 天 牛 科 中 。 虽 然 本 实验 成 功 扩 增 出 了 63 种 天 
牛 18S rDNA VA 和 V7 区 以 及 28S rDNA D2 和 D3 
区 的 序列 ,然而 在 亚 科 水 平 上 仍 需要 补充 做 更 广泛 
的 研究 ,为 解决 天 牛 总 科 高 阶 元 进化 关系 问题 提供 
依据 。 
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